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1 INTRODUCTION 



The discovery of cosm ic acceleration, (|Riess et al.l 1 19981 : 
iPerlmutter et af]|l999t ). and its subsequent explanation as 
dark energy has resulted in the proposal of many alternate 
models to the current cosmological paradigm A Cold Dark 
Matter (ACDM) as the nature of dark energy remains elu- 
sive. There are three main alternatives to dark energy: either 
General Relativity is not the correct description of gravity 
on large scales, introduce a higher dimensional theory, or 
the metric that describes t he universe is not required to be 
maximally symmetric (see iDurrer fc Maar tens '*20 0^, for a 
review). The Fractal Bubble (FB) model (|Wihshirc 2007al lbh 
falls into the last category, as dark energy is replaced by 
'quasi-local gravitational energy' arising from changes in the 
curvature of spac etime due to an in homogeneous distribu- 
tion of matter. In lLeith et al.l l|2008h . it is claimed that the 
FB model achieves a concorda nt set of cosmologic al parame- 
ters using supernova data from lRiess et al.l (120071'). the angu- 
lar sc ale of the sound horizon from WMAPl (iBennett et al.l 
l2003l ) and the com oving spatial separatio n of the correlation 
function in SDSS (|Eisenstein et all 120051 ). 



ABSTRACT 

The Fractal Bubble model has been proposed as a viable cosmology that does not 
require dark energy to account for cosmic acceleration, but rather attributes its ob- 
servational signature to the formation of structure. In this paper it is demonstrated 
that, in contrast to previous findings, this model is not a good fit to cosmological 
supernovae data; there is significant tension in the best fit parameters obtained from 
different samples, whereas ACDM is able to fit all datasets consistently. Furthermore, 
the concordance between galaxy clustering scales and data from the cosmic microwave 
background is not achieved with the most recent supernova compilations. The valid- 
ity of the FB formalism as a sound cosmological model is further challenged as it is 
shown that previous studies of this model achieve concordance by requiring a value 
for the present day Hubble constant that is derived from supernovae data containing 
an arbitrary distance normalisation. 
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Although the vast majority of inhomogeneous cos- 
mological models have been shown to be inconsistent 

with observa tions or rely on a p a rticular cho ice of 

coordinates lllshibashi fc Waldl l2006f : IZiblin et all l2008l : 
ICaldwell fc Stebbins II2008I ). the FB model remains unchal- 
lenged. In this paper, we examine the ability of the FB 
model to produce a concordant set of cosmolog ical param- 
eters. In particular, we address the claim of iLeith et al.l 
(|2008h that their two parameter fits for the distance mod- 
ulus given by the FB model us i ng th e supernova type la 
(SNe la) sample of [R iess et all (|2007l ) supports a concor- 
dance of observational evidence for the FB model. Addi- 
tionally, we demonstrate a number of discrepancies includ- 
ing the i nability of t h e FB model to consistently describe 
both the iRiess et al.' ^200?) sample and the Union super- 
nov a compilation (Ko walski ct al. 2008.) or the Constitution 
set (|Hicken et al.ll20o'^) despit e none of these SNe la samples 
bei ng in tensio n. Leith et al.l l|2008l ) also omit to calibrate 
the[R iess et ahl (|2007l ) SNe la data to a distance scale when 
fitting for the Hubble constant, Hq. 
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2 BACKGROUND 

The FB model (|Wihshirdl2007al lbl) is a two scale model with 
a local regime for over or underdense regions and another 
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Figure 1. Top: Marginalised post erior distrib utions for the FB 
model using the Union compilation llKowalski et al. 2008) in solid 
lines, the Gold SNe la sample from Riess ct al. (2007) in dot- 
ted lines, and the Constitution set (Hicken ct al. 2009) in dashed 
lines. The colours correspond to the confidence limits, red is Icr, 
green is 2ct and blue is Scr. Bottom: As above except using ACDM. 



for the volu me average, w hich utilises the Buchert averag- 
ing scheme (jBuchertI l2000l ) . We use an overbar to indicate 
volume averaged quantities as opposed to 'dressed' param- 
eters measured by observers in galaxies, termed 'walls' in 
the FB model. The volume averaged scale factor, a, and the 
void fraction, /„, defined as the total volume in void regions, 
evolve according to the Buchert equations as follows; 
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where a is the normalisation constant of the cur- 
vature energy. Equation 1 is equivalent to writing: 
— QqH^ — QkH^ = QrriH^, sincc the volume averaged 
Hubble constant is defined as H — d/a and the normalised 
volume averaged energy densities are given by: 
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3 d^m 

Information about the distribution of matter is encoded 
in the void fraction, /„, which is directly proportional to 
the dressed normalise d matter density , f2m,o, as 0^,0 ~ 
1(1 - /«,o)(2 + /„,o) l| Wiltshire! [2007^ ). Although a back- 
reaction term occurs in the above equations, this is not the 
principal mechanism that creates cosmic acceleration, but 



rather this is achieved through 'quasi-local' gravitational en- 
ergy quantified by the curvature energy term. Such gravita- 
tional energy can not be localised in the stress energy tensor 
and in the FB model gives rise to an apparent cosmic ac- 
celeration as measurements are distorted when signals pass 
through regions of varying curvature and it cannot be as- 
sumed that the clocks of all observers are synchronous inde- 
pendent of location. In the absence of homogeneity, volume 
averaged quantities must be transformed back to what ob- 
servers in overdensities like ourselves would measure. The 
FB model, chooses to foliate spacetime such that there is 
a lapse function, defined as 7(t) = dt/dr, where r is the 
proper time of an observer in a wall or a void and t is the 
volume averaged or cosmic time, which quantifies the dif- 
ference between clocks of two different observers depending 
on their location. If the lapse function in walls is signifi- 
cant, then the FB model predicts that a spurious cosmic 
acceleration will be detected by failing to account for an 
inhomogeneous matter distribution through the assumption 
of a Friedmann-Lemaitre-Robertson- Walker (FLRW) back- 
ground and a recalibration of observations for the difference 
in clocks between walls and the volume average is necessary. 
This difference may be as l arge as 7 1.38 t o achieve the 
level of concordance seen in iLeith et all l|2008l ). 

Unlike previous approaches with inhomogeneous cos- 
mologies, much of the physics behind the FB model is es- 
sentially alien to the standard model of moder n cosmol- 
ogy. An extension to the equivalence principle l|Wiltshirel 
l2008h is required to explain the anomalously large value 
of 7 necessary to pr oduc e the best fit par ameters quoted 
in iLeith et ail (120081 ) and IWihshiri (|2007a). Furthermore, 
neither the Newtonian limit or Birkhoff 's theorem is relevant 
in this cosmology. Leaving aside these conceptual issues, we 
have instead focused on establishing the compatibility of 
currently available observational data with the FB model. 



3 TENSION WITH SUPERNOVA DATA 
3.1 Constraints from new data sets 



The Union supernova compilation (|Kowalski et al.ll2008l ) is 
primarily composed of SNe la catalogs from literature that 
have been reanalysed in an uniform manner to reduce sys- 
tematics, with some further cuts imposed to exclude data 
of insufficient quality. The compilation contains a total o f 
307 SNe la, of which 27 are drawn from iRiess et al.l (|2007D 
after the final cut and a further eight are derived from new 
low-redshift observations. Its constituent catalogues have 
also been shown to be consistent with one another und er 
this new analysis [see Figure 9 of iKowalskT et al.1 (|2008l )]. 
More recently, the entire Union compilation has been amal- 
gamated with an additional sample of 90 new CfA3 SNe 
la into the Constitution set, incr easing the amount of low 
redshift data by a factor of 2.6 (|Hicken et al.|[2009l ). Like 
the Union compilations, the Constitution set is uniformly 
reduced through the use of a single light curve fitter for all 
SNe la. Clearly, both the U nion compilation and the Gold 
sample of lRiess et al.l (|2007l ) are entirely consistent with the 
Constitution set. 

We test the ability of the FB model to constrain cos- 
mology with the new Union and Constitution SNe la data 
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Figure 2. Percentage of SNe la from Gold (solid), Union (dotted) 
and Constitution (dashed) samples falling into redshift bins of 
width Az = 0.2. There are 397 SNe la in the Constitution set, 
307 SNe la in the Union com pilation and 182 SNe la in the Gold 
data set of lRiess et al ] l|2007h . 



sets, as well as the Gold sample f rom Riess et al. (^2007^ 



This last sample was analysed by iLeith et all (|2008l ). but 



this study suffered from a number of problems as, detailed 
in Section 3.2, which we rectify in this work. The best fit 
cosmological parameters are found from calculating the re- 
duced chi squared, x^, based on the distance modulus given 
by /X = 51ogio(dL)+25 for ACDM and FB model to aU three 
samples of SNe la. The luminosity distance d^, is defined as: 

(l + -)«o/;°^ (ACDM), 
i.„n ^ _ /„)i/3 j;o (FB). 



dL = 



i«o(l + ^)(l 



7(l-/^.)l/^^^ 

There are two free parameters of interest in the FB model 
that are used to calculate the distance modulus; the cos- 
mology is fully constrained after specifying the initial void 
fraction the dressed Hubble constant today. Ho, and 
the initial ratio of the expansion rates as measured in a wall 
and a void, hri\ the latte r remains fixed for all fits. We have 
used the same priors as iLeith et al.l (|2008l ) on f„i and h^i, 
namely lO"'^ < f„i < 10"^ (and hence 0.01 < f2™,o < 0.95) 
and hri = 0.99999 at last scattering, with 50 5S Ho ^ 75 
kms~^Mpc^^. The ACDM priors are the same on Ho and 
cover a similar range on Q.mfi'- ^ fim.o 1- The value of 
Ho is treated as a nuisance parameter, which we marginalise 
over, since it cannot be generally assumed that the absolute 
magnitudes of the SNe la are known. This results is a single 
parameter fit for r^m.o, which is directly related to fv 

The statistic is converted to a likelihood, £, via the 
relationship, C oc exp(— x^/2). Bayesian inference then gives 
the resultant posterior distributions as presented in Fig- 
ure [1] for both cosmological models, with the best fit val- 
ues shown in Table [T] For ACDM, all three distributions 
give a consistent set of cosmological parameters that are 
in agreement, with the la confidence regions significantly 
overlapping. This does not occur for the FB model; the la 
limits derived for each SNe la sample do not all coincide. 
We must consider the 2a confidence region before we can 
find a value of f2,„,o that agrees with between the Union 
co mpilation o r the Constitution set and the Gold sample 
ofE less et ahl (|2007|). In fact, t he best fit value for the Gold 
sample from lRiess et al.l l|2007l ) is ruled out at So" when per- 
forming the same analysis for the Union compilation and is 
beyond the 3a limit for the Constitution set. (See Table [1] 



for the complete set of best fit parameters and la errors). 
However, there is good agreement between the best fit FB 
parameters for Constitution set and the Union compilation, 
which is reassuring; the reason for this discordance lies with 
the model rather than an anomaly in the SNe la sample. If 
this model is to be acknowledged as a viable alternative to 
ACDM, then it seems that we must also accept that there is 
no consistent observational evidence from SNe la for the FB 
mo del or that t h ese tw o SNe la samples are in tension with 
the Riess et al.l l|2007h data set. However, iKowalski et al.l 
(|2008l ) found a high degree of consistency between the sam- 
ples used in the construction of t he Union compilat i on, albe it 
with mild tension betw e en th e IKrisciunas et all (|2004al lbl. 
I2OOII) and iHamuv et al.l l| 19961 ) samples and the other SNe 
la data when comparing dn/dz. Crucially, the more recent 
and populous samples break the concordance of cosmolog- 
ical tests for the FB model; clearly the best fit values for 
f2m,o as shown in Table [T] do not overlap with those ob- 
tained from galaxy clustering s tatistics and the an gular scale 
of the sound horizon found bv iLeith" et all (|2008l ). which re- 
quire .27 < Om.o S , 0.37. Furthermore, although the Bayes 
fac tor (lTrottall2qod) b etween the FB model and ACDM with 



the lRiess et al. I 2007[) Gold data set ma rginally favours the 
former [InB = 0.27, iLeith et all l|2008l )]. the FB model is 
weakly disfavoured when using the Union compilation or 
the Constitution set (InB = -1.38, -1.469, respectively) un- 
der the same classification scheme on the strength of evi- 
dence. As the treatment of systematics improve and sample 
sizes grow, ACDM is increasing being favoured by the data 
while the best fit parameters are moving away from those 
required by the FB model for concordance. The deteriora- 
tion of the evidence towards the FB model in terms of the 
Bayes factor may be attributed to the width of the posterior 
in flm,o, which remains significantly larger than that of the 
corresponding data set under ACDM. 



3.2 Distance normalisation? 



Figure 2 of iLeith et all l|2008l ) depicts a concordance dia- 
gram, represented in the fi^.o-Ho plane, for the FB model, 
w ith confidence lim its from SNe la using the Gold data 
of lRiess et all (|20o3), the proper distance to the sound hori- 
zon seen in SDSS, and the angular scale of the sound hori- 
zon at decoupling. However, a number of issues regard- 
ing the analysis of these observational constraints remain 
unaddressed by the authors, which, when corrected, fal- 
sify the concordance. The most serious of these is that the 
SNe la data used to produce this diagram contains an ar- 
bitrary value of Ho without the appropriate correction to 
the distance moduli. Thus, it is misleadin g to pr esent con- 
fidence limits in Figure 2 of iLeith et aU (|2008D that de- 
pend on Ho, since these do not repre sent its true value. 
Instead, the concordance diagram in iLeith et al.l (|2008l ) 
should be shown with either Ho as a nuisance parameter, 
whic h is marg i nalise d or with the recommended calibration 
of , Riess et al.ll|2007h . 

For the Gold sample, iRiess et "ahl (|2007l ) suggest us- 
ing a systematic subtraction of 0.32 magj to use the same 



^ Magnitude calibration and data can be obtained online at 
br aeburn . pha . j hu . edu/ ~ariess/R06/ sn_ sample 
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Figure 3. Le/t: Residuals from an empty FLRW universe with 
same value of Ho as the FB model with best fit parameters from 
the Union compilation ( solid), Constitution set (dot-dashed) and 
Gold iRiess et al.l l|2007l ) sample (dashed) . The diffe rence in the 
distan ce moduli between the best fit values found bv lLeith et alj 
l|2008l ) and those from the more recent compilations are shown 
in light grey (Union) and dark grey (Constitution) Right: 
Range of residuals from an empty FLRW universe with same 
value of Ho for FB model (solid grey) and ACDM (hatched 
black) for 0.1 < f!m,o < 0.5. Any residual from an empty 
universe using either the FB model or ACDM, with the above 
parameters, will fall into the corresponding shaded area. The 
curves show the difference between the distance moduli of the 
FB model and ACDM for the same value of Hq for m , o — 
0.30 (solid), 0.35 (dashed), 0.40 (dot - dashed), 0.45 (dotted). 
There is a preferred range of parameters, namely 0.35 ^ Q-m. ,0 < 
0.4 for which the difference between the two models is minimised. 



Cepheid scale as iRiess et al.l (|2005l ) if the value of Ho is of 
interest and this is unnecessary for fits to dynamic quanti- 
ties only. An arbitrarily chosen value of Ho is inserted into 
the data, since the absolute magnitude of a SNe la is not 
known until an appropriate distance scale, such as that ob- 
tained from Cepheid luminosities, is applied. This correc- 
tio n to the distance m oduli has not been taken into account 
bv lLeith et al. (|2008|), who claim to follow the Cepheid cal- 
ibration of ISandage et al] (|2006| ) but no corrections have 
bee n applied a t all. Indeed, when performing our fits to 
the iRiess et al.l (|2007l ) Gold sample, we are ab le to repro- 
duce t he same best fit value of Hq as quoted in iLeith et al.l 
(|2008l ) (within 0.1 error) without accounting for the arbi- 
trary distance normalisation at all. The authors, however , 
have interpreted the value of Ho in the iRiess et al.l (|2007l ) 
sample as a physical parameter t hat constrai n s cosm ology, 
rather than as a value ch osen by Riess et al.l ([20o3). Thus 
the fit to Hq in Figure 2 of iLeith et al.l ( 20081 ) is meaningless 
without calibrating the data. It is plausible th at the arbi- 
trary normalisation chosen bv lRiess et al.l (|2007l ) happens to 



coincide with the Sandage distance scale, but no justificatio n 
is given for the choice of calibration in iLeith et al.l l|2008h . 
Int erestingly, using t he Cepheid calibration recommended 
bv I Riess et al.l l|2007l ) gives a best fit value of Ho = 71.6 
kms~^Mpc~^ with Qm,o unchanged, since changing Ho only 
shifts the scale of the Hubble diagram. This no longer pro- 
duces a conco rdant set of cosmo logical parameters; in fact 
it is stated in iLeith et al.l (|2008h that for any value of Ho 
greater than 70 kms~^Mpc~^, the observational constraints 
from the proper distance to the sound horizon and its angu- 
lar scale at decoupling would no longer agree. 



4 A COINCIDENCE PROBLEM FOR THE FB 
MODEL 

To understand why the cosmology of the FB model changes 
dramatically for the Union compilation and the Constitu- 
tion set while ACDM remains consistent, we have analysed 
the redshift distributions of these two SNe la cata l ogues . 
Figure [5] shows that the Gold sample of iRiess et al] (|2007l ) 
contains a slightly higher percentage of high redshift SNe la 
than in the other samples but there are significantly more 
low redshift samples in the latter, particularly in the 0.0 
^ 2 <0.2 range for the Constitution set and 0.35 < z < 0.4 
range for the Union compilation. Furthermore, it should be 
unsurprising that the FB model is able to fit both the Con- 
stitution set and the Union compilation within la error; 
Figure [2] also shows that the redshift distribution of these 
two samples are very similar, especially in the high- 2: regime 
where the Constitution set is dominated by data taken from 
the Union compilation. 

The reason behind the shift in parameter space is made 
apparent by considering the range over which the newer 
compilations contain more data. This is coincident with 
the region in which the three different best fits for the 
FB model diverge the most (Figure O left). The residu- 
als have been taken relative to an empty universe with 
^m,o ~ 0, f^A.o = 0, flk.o = 1, and the same Ho to remove 
the uncertainty in its value. For the Union and Constitution 
sets, the FB model requires a much greater amount of ac- 
celeration to fit the data and indeed the period of apparent 
acceleration occurs at 2 ~ 1.65 and 2 ~ 1.85 respectively. 
As the cosmological constraints become progressively tighter 
in the more recent SNe la samples, the FB model is increas- 
ingly pushed to higher values of curvature energy. Indeed, 
the best fit parameters from the Union and Constitution 
sets imply that the universe is almost entirely composed of 
voids with an almost negligible amount of collapsed struc- 
ture. As remarked in the previous section, the FB model has 
a flatter, wider posterior than ACDM, which implies that a 
greater range of the parameter space can fit the data. This 
can be attributed to the limited range of behaviour that the 
FB model can exhibit in terms of changes to the distance 
modulus when /„i is altered. This is demonstrated in the dif- 
ference between the range that is spanned by the residuals 
for the FB model (solid grey) and ACDM (hatched black) 
for 0.1 < Qm,o < 0.5 (Figure O right). Although the prior 
on /„i spans ~ 3 orders of magnitude, the residuals for the 
FB model have a substantially truncated regime of accel- 
eration in comparison to ACDM. This partly accounts for 
the discordance seen in the confidence limits for this model 
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Table 1. Best fit parameters for ACDM and Fractal Bubble 
model to SNe la data. Distance normalisation is arbitrary for 
the Union compilation and the Constitution set. Best fit values 
of Om,o S'lMi Ho are subject to an additional error of ±5 X 10^^ 
and ±0.1 respectively due to grid spacing. 



Model 



SNe la sample 



Ho (km/s/Mpc) 



FB 

ACDM 
FB 

ACDM 
FB 

ACDM 
FB 

ACDM 



Riess et al. (2007)^ 
Riess et al. (2007)^ 
Riess et al. (2007') ^ 
Riess et al. (2007*) ^ 
Kowalski et al. (2008*) 
Kowalski et al. (2008) 
Hicken et al. (2009) 
Hicken et al. (2009) 



0.333 
0.34j 
0.33 
0.34 
0.08 
0.29 
0.02j 
0.29^ 



+0.07 
-0.11 
+0.045 
0.035 
+0.08 
0.08 
+0.03 
-0.025 



61.8 

62.8 
71.6 
72.8 



^Without calibration 
^With calibration 



that is not replicated in ACDM. The best fit parameters of 
the FB model are extremely sensitive to small changes in 
the SNe la data as it needs to compensate for these by a 
large variation in /„i when fitted to an another redshift dis- 
tribution with a different amount of error on each SNe la. 
In addition, there is a special set of values for fvi which will 
mimic ACDM parameters well, that is the dressed matter 
density, f2m,o, in the FB model is a similar value to Slm.o 
derived from ACDM. 

In the right panel of Figure O we have also considered 
the difference in distance moduli between ACDM and FB 
models (each with the same value of Ho). This is minimised 
when 0.35 < f2m,o ^ 0.4 for both models. While it is strictly 
not necessary for the FB model to predict the same value of 
i^m.o as ACDM, since its value is not directly observable but 
in ferred f rom theory, it does so because the SNe la sample 
of lRiess e t al. (2007) gives a best fit value of Q.m,o ~ 0.34 for 
ACDM and this happens to fall within this parameter range. 
Since the best fit parameters for both the Union and Con- 
stitution set are f2m,a = 0.29 and Q.a,o ~ 0.71 for ACDM, 
fitting the FB model to this dataset cannot produce a sim- 
ilar value for flm,o- In addition, changes in the data that 
are insignificant for CDM cosmologies will have sufficient 
leverage to skew the posterior of FB models as a result of 
the similarity of the solutions when fvi is varied. This be- 
haviour can__beattributed to the existence of a tracker so- 
lution l|Wiltshirel |2007b) such that all solutions converge at 
high redshift, which inhibits the range of possible behaviour 
that the FB model can exhibit. 



5 CONCLUSIONS 

One of the more problematic features of the FB model 
is its failure to provide any predictions on cosmological 
parameters that are directly observable. Although it pro- 
duces a val ue of Ho that is con sistent with the EST 
Key Project (|Freedman et al. | 2001 ) and obse rvations from 
WMAP when fitted to the lRiess et al.l (120071 ) SNe la sam- 
ple with the appropriate Cepheid calibration, this value is 
not concordant with that obtained from the proper distance 



to th e sound horizon observed in SDSS (jEisenstein et al.] 
l2005h or the angular scale of the sound horizon observed 
in WMAPl (|Bennett et al.l 120031 '). Furthermore, the pre- 
dictions that are made require further investigation of the 
model at a fundamental level; the constrains from differing 
SNe la sample are discordant . Regardless of the calibration 
chosen for the Gold sample of lRiess et al.l (|2007h . there is no 
value of Ho that can produce a consistent cosmology within 
the FB framework under the currently available SNe la data. 

The most appealing feature of the FB model is to offer 
a mechanism for replacing dark energy with an inhomoge- 
neous matter distribution and yet it lacks any formalism 
by which structure formation can occur. Despite the void 
fraction being a fundamental parameter of the FB model , it 
ha s not been made tr ansparent in either [Wiltshire! (|2007a| [bh 
or iLeith et al.l l|2008l ) how this is manifested in measure- 
ments of clustering statistics or the growth of structure. In- 
de ed, the observation al evidence for the FB model presented 
bv lLeith" et al.l l|2008l ) are all derived from geometrical tests 
and it is not clear what behaviour it would exhibit under 
more dynamically oriented probes such as large scale struc- 
ture surveys or constraints from weak lensing. Although it 
is claimed in lLeith et al.l (|2008| ) that the 'concordance' val- 
ues of ilm,Q derived from the FB model are more compati- 
ble wit h those obtained fr om X-ray measurements of cluster 
counts (|Yepes et al.ll200"7l ). it is in fact impossible to state 
what the FB model would predict without a mechanism for 
structure formation, since such estimates of ilm,o are model 
dependent. Until the FB model can provide more plausible 
observational evidence, it is difficult to envisage this model 
as a serious competitor to ACDM when so many questions 
remain. 
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